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^\ ' Hybrid Chaplygin gas model is put forward, in which the gases play the role of dark energy. For 

this model the coincidence problem is greatly alleviated. The effective equation of state of the dark 
energy may cross the phantom divide w = — 1. Furthermore, the crossing behaviour is decoupled 
from any gravity theories. In the present model, w < — 1 is only a transient behaviour. There is 
(— | a de Sitter attractor in the future infinity. Hence, the big rip singularity, which often afflicts the 

models with matter whose effective equation of state less than —1, is naturally disappear. There 
I—} ' exist stable scaling solutions, both at the early universe and the late universe. We discuss the 

| perturbation growth of this model. We find that the index is consistent with observations. 

PACS numbers: 95.36.+X 
Ck I. INTRODUCTION 

2 : The existence of dark energy is one of the most significant ecological discoveries over the last century Q |. 

Although fundamental for our understanding of the Universe, its nature remains a completely open question nowadays. 
, ^ , ■ For recent studies of dark energy, see review article Q . 

Recently the so-called Chaplygin gas, also dubbed quartessence, was suggested as a candidate of a unified model 
of dark energy and dark matter The Chaplygin gas is characterized by an exotic equation of state (EOS) 

> 



Pch = -A/ p ch , (1) 

where A is a positive constant. The above equation of state leads to a density evolution in the form 



^ = y/A + B(l + z)*, (2) 

where A = Aj p\, B is an integration constant, z denotes the redshift, po represents the present critical density. The 
L" ■ attractive feature of the model is that it naturally unifies dark energy and dark matter. The reason is that, from |2|). 
. £^ \ the Chaplygin gas behaves as dust-like matter at early stage and as a cosmological constant at later stage. Some 
possible motivations for this model from the field theory points of view are investigated in [1, 0| . The Chaplygin gas 
emerges as an effective fluid associated with ci-branes 0, Q and can be also obtained from the Born-Infeld action [9[ . 

The Chaplygin gas model has been thoroughly investigated for its impact on the 0th order cosmology, i.e., the cosmic 
expansion history (quantified by the Hubble parameter H [z]) and corresponding spacetime-geometric observables. An 
interesting range of models was found to be consistent with SN la data [T5L Till [l2l [l3| , CMB peak locations 
and gas mass fractions in clusters of galaxies [l5j . There seems to be, however, a flaw in unified dark matter (UDM) 
models that manifests itself only on small (galactic) scales and that has not been revealed by the studies involving 
only background tests. In [16], it is found that generalized Chaplygin gas (GCG) model produces oscillations or 
exponential blowup of the matter power spectrum inconsistent with observations. In fact, from this analysis, 99.999 
% of previously allowed parameter of GCG model has been excluded (see, however, fl7l. Il8l. Il9j ) . 

Hence we may turn to a model with Chaplygin gas and dark matter. It has been pointed out that Chaplygin 
gas model can be described by a quintessence filed with well-connected potential Q. Therefore, a model with 
Chaplygin gas and dark matter is essentially a special quintessence model, in which the EOS of dark energy (Chaplygin 
gas) w (defined as the ratio of pressure to energy density) satisfies < w < — 1. 
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A merit of the Chaplygin gas model, in which Chaplygin gas only plays the role of dark energy, is that the 
coincidence problem is greatly alleviated. The coincidence problem in ACDM model says that A keeps a constant 
while the density of CDM evolves as (1 + z) 3 in the history of the universe, then why do they approximately equal 
each other at "our era"? We see from ^ that density of Chaplygin gas evolves as the same of dark matter and 
only at late time it evolves as cosmological constant. Therefore, the coincidence problem may be alleviated at some 
degree. We shall discuss this possibility. 

The cosmological constant is the far simple candidate for dark energy. Though previous observations are consistent 
with the cosmological constant , they leave enough space for a dynamical dark energy [20|, |2l|, |22j, |23j . Following 
the more accurate data, the implications for dynamical dark energy become clear: the recent analysis of the type la 
supernovae data indicate that the time varying dark energy gives a better fit than a cosmological constant, and in 
particular, the equation of state parameter w crosses —1 at some low redshift region from above to below 23, 25 . l26l. 

I23, Is^, 12^, l^iL laS, laS, 1^3, laa laE> la^, la^, E3, E3, 1^^, S3, Eg, S3, 1^3, IHSi IHM, IBS, IHa IHa - it 



seems difficult to realize this transition in context of Chaplygin gas dark energy model. However, if we consider EOS 
(1) carefully and rigorously, we shall find the other solution of the continuity equation different from (2), 

Pch 



-^C + D{l + zf, (3) 
Pa 



which satisfies the EOS, 



C Pl 



2 



PCh = -^, (4) 
pch 

where C and D are constants. A noticeable property of this case is that the energy density of Chaplygin gas is 
negative. But this is not anything completely new. If quantum effect is considered, the energy density of a field can 
be naturally lower than zero [571 ]. Also in the phantom dark energy model which is extensively studied in context of 
cosmology [58[ , the density of dark energy can be less than zero if observed by observers other than the homogeneous 
observers, since it violates the weak energy condition. We call this Chaplygin gas with negative energy density 
"type II" Chaplygin gas , and correspondingly, the original Chaplygin gas with positive density is dubbed "type I" 
Chaplygin gas in the present article. Therefore, we find there are two branches of mathematical solution satisfying 
the EOS of Chaplygin gas and continuity equation. In fact, B in (2), as an integration constant, can also be negative. 
However, a difficult arises when B < 0: the universe bounces at some finite redshift. To avoid to plague the success 
of nucleosynthesis, < 6 x 10~ 28 , which is an unnatural fine-tuning. In this paper we constrain ourself in the 

case of B > 0. 

A kind of interacting Chaplygin gas model in which the Chaplygin gas plays the role of dark energy and interacts 
with cold dark matter particles has been investigated in [59[ . In this model the effective equation of state of Chaplygin 
gas may cross the phantom divide. 

We shall present a hybrid model composed by both of the two types of Chaplygin gas, in which the EOS of 
dark energy can cross the phantom divide, which is decoupled from any gravity theory. Assuming standard general 
relativity, we prove there exist scaling solutions, both in the early universe and the late time universe: the early 
universe is attracted by a dust tractor, that is, the cosmic fluids enter a dust-like phase; and the late time universe is 
attracted by a de Sitter tractor , that is, it evolves into a de Sitter phase. 

Often, the bi g ri p singularity is narrated as follow: at finite cosmic time, both the scale factor and energy density 
become infinite [6fj|. This is a conclusion in a special chat, although in the most used chat. As we learned from 
Schwarzschild solution, an ordinary point can be a singularity in disguise for a special chat. We shall present some 
discussions on this issue and show the big rip singularity is a true singularity. The big rip singularity usually emerges 
at the models with phantom like dark energy, because the density of a phantom field will increase with time. But in 
our present model w < — 1 is only a transient behavior, and the final state of the universe is always a de Sitter, for 
the whole parameter space. 

To construct a model simulating the accelerated expansion is not very difficult. That is the reason why we have 
so many different models. Recently, some suggestions are presented that perturbation (fluctuation) growth function 
5{z) = 8p m / Pm of the linear matter density contrast as a function of redshift z can be an effective probe to explore 
the dark energy models 0, M, M, The related works can be traced back 

to SSSlzI. 

There is an approximate relation between the growth function and the partition of dust matter [6l|, [62|, HH, [H, l65l . 

M M HI M m M M M M , 
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where Q, m is the density partition of dust matter, a denotes the scale factor, and 7 is the growth index. This 
relation is a perfect a ppr oximation at high redshift region. Also, it can be used in low redshift region, see for example 
ill, EI EI IB! EI EI M, El EI H IzE El IZI III . The theoretical value of 7 for ACDM model is 6/ 1 1 [zl M, [Hi III ■ 
We shall investigate the perturbation growth of the present model. 

We present our model in details in the next section. In section III, we study the perturbation growth in the hybrid 
Chaplygin gas model. Our conclusions and discussions appear in the last section. 



II. THE HYBRID MODEL 



In this section we shall investigate three properties of the hybrid Chaplygin gas model. In the first subsection 
we display that the coincidence problem is greatly alleviated. In subsection B, we study the crossing —1 behavior of 
the EOS of the dark energy. And we explore the the evolution of this model by the dynamical-system analyses in 
subsection C. 



A. Coincidence problem 



We consider a model in which the type I Chaplygin gas and type II Chaplygin gas play the role of dark energy 
together. The density of dark energy reads 

Pde = Pl+P2- (6) 

The continuity equation for which reads, 

dp de = Z(p de +p de )d\\i(l + z). (7) 

Here pi denotes the density of type I Chaplygin gas and P2 represents the density of type II Chaplygin gas , 

Pi 
Po 

P 

The energy density of dust matter redshifts as 



VaTW+W, (8) 

Po 



w — = -y/C + D(l + z)e. (9) 
Po 

p^ = n mQ (i + z) 3 . (10) 

Po 

Hence, if dark energy is just cosmological constant, we suffer from a coincidence problem. However, if Chaplygin 
gas plays the role of dark energy, this problem will be alleviated because at the early universe the Chaplygin gas also 
behaves as dust. 

We see from figure Q] that if a cosmological constant plays the role of dark energy, the density of matter is about 
20 orders larger than the dark energy at z = 1100: they must have been fine-tuned at that time. However, if the 
Chaplygin gas serves as dark energy, the ratio of dark matter and dark energy keeps at order 1, which relieves the 
coincidence problem, for all 3 cases of the universe. 



B. Crossing —1 

From the continuity equation ([7]), we arrive at 

Pde , 1 d\np de 
Pde 3dln(l + z) 

which means that in an expanding universe if p de decreases and then increases with respect to the redshift, or increases 
and then decreases, then we conclude that EOS of dark energy crosses phantom divide. Here, 

^ = -3[Bp^(l + zf + Dp- 2 \\ + zf]. (12) 

We see that if we carefully tune B and D, the EOS of dark energy may cross —1. In figure [2l we show some concrete 
examples for which Chaplygin gas crosses the phantom divide at about z = 0.2, where a = p de /po- Also, in figure[3] 
we show the phase portrait of a vs. j3 = p m / P01 where p m represents the density of dust matter. A note is that any 
gravity theory is not interposed up to now. Our results only depend on the continuity equation. 
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FIG. 1: A double logarithmic plot: the densities of dark energy as functions of redshift. The cosmological constant inhabits on 
the horizontal line, the dust matter dwells on the thick solid curve, and other three curves represent the Chaplygin gas dark 
energy. The long dashed, short dashed, and thin solid curves denote the dark energy in a positive curvature, negative curvature 
and flat universe, respectively. 
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FIG. 2: a as a function of z. On all of the 5 orbits B = 
the below A = 1.104, C = 0.048; A = 1.059, C = 0.039; 



0.3, D = 0.15, n m = 
A = 1.01, C = 0.03: 



0.28 with different A and C. From the above to 
A = 0.9106, C = 0.015; A = 0.8267, C = 0.006 



C. Dynamical analysis 

Often, a model containing phantom-like matter is afflicted by big rip problem, that is, at some finite time, the scale 
factor and density is divergent. However, this description of singularity depends on FRW coordinates, which may be 
an ordinary point described by other coordinates. We check scalarpolynomial curvature, which is a good lesson we 
learned from Schwarzschild solution. Here we take the example in [79(. Consider a universe is dominated by dust and 
then by constant w phantom. The transition occurs at t pm . At 



wt 



pm 



1 



W 



(13) 



the scale factor and density become divergent. The simplest scalar polynomial is Ricci scalar i?, 



R 



16 4 

OL pm L m 



(1 + w) 



-w + (1 + w)- 



(14) 
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FIG. 3: Phase portrait for a vs. (3. On all of the 5 orbits B = 0.3, D = 0.15, £l m = 0.28 with different A and C. From 
the above to the below A = 1.104, C = 0.048; A = 1.059, C = 0.039; A = 1.01, C = 0.03; A = 0.9106, C = 0.015; A = 
0.8267, C = 0.006 

Because w < — 1, the term in first square bracket always larger than zero, and term in the second bracket is divergent 
when t — > w x *™ . Therefore, the scalar curvature is infinite when it goes to the big rip point. Also, because all 
comoving observers move along timelike geodesies, and the cosmic time is their proper time, this singularity is also an 
incomplete-geodesic singularity. Hence, big rip singularity is a true singularity. Through researches on the dynamical 
properties of the universe with Chaplygin gas and dark matter, we shall show that there is no future singularity in 
this model, though the dark energy behaves as phantom in some stage. 

Now we start to study the dynamical evolution of the universe in frame of standard general relativity, for which we 
introduce Friedmann equation in a spatially flat FRW universe, which is implied either by theoretical side (inflation 
in the early universe) ,or observation side (CMB fluctuations [231] ^ . 

H 2 = 3^2 (P1+ P2+ Pm) (15) 

where, as usual, H denotes the Hubble parameter, and /j, stands for the reduced Planck mass. For convenience we 
first define the following new dimensionless variables, 

x 

y 

u 

h 
h 

The dynamics of the universe can be described by the following dynamical system with these new dimensionless 




(16) 
(17) 
(18) 
(19) 
(20) 
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variables, 

x> = -*- x -\ x *- x -Ht) + \xP, (21) 

y' = -\y-\y 2 -y- 2 li) + lyP, (22) 

u = -\u+\uP, (23) 

l'i = \hP, (24) 

1' 2 = \hP, (25) 

where 

P = x 2 ~x- 2 li-(y 2 -y- 2 4) + u 2 , (26) 



and a prime stands for derivation with respect to s = — ln(l + z). Note that the 4 equations (|2ip 
of this system are not independent. By using the Friedmann constraint, which can be derived from the Friedmann 
equation, 

x 2 ~y 2 + u 2 = l, (27) 

the number of the independent equations can be reduced to 4. The critical points of this system satisfying x' = y' = 
I' = b' = appearing at 

h=h= 0, (28) 
x 2 -y 2 + u 2 = 1, (29) 

and 



u = 0, (30) 
x = h, (31) 
y = h- (32) 

The first set of critical points (|28p . (|29|) dwells at the early universe, since H — > oo. Also this set satisfies the constraint 
equation automatically. One sees that it is fairly ample, which only needs (x c ,y c ,u c ) inhabits on the surface (|29l ). 
where c label the critical point. The reason roots in the fact that all of the three components, type I Chaplygin gas, 
type II Chaplygin gas and matter are dust-like in the early universe, they can evolve into a scaling solution with rather 
arbitrary proportion of components along inverse time direction. We call this set of critical points dust attractor. The 
second set of critical points resides at the late time universe, because u = 0, which means matter has been infinitely 
diluted. By using constraint equation (|27[) . we further derive x c — y c = l, which resides on a hyperbola. Because this 
set of critical points ensure that l\ = I2 = constant, which means H = constant, we call it de Sitter attractor. The 
previous models in which the EOS of dark energy crosses —1 often suffer from future or past difficulties, such as big 
rip disaster or to plague the structure formation theory because of a too stiff EOS. Hence, most of them only can be 
used to describe the evolution of the universe at some low redshift. By striking contrast, our model are free of such 
difficulties from CMB decoupling to the future infinity. Global fittings or simulations of structure formation operate 
routinely in frame of the present model. 

Though in above context we call the singularity "attractor" , it is only an intuitive conclusion. To obtain a mathe- 
matically strict result of the stability at the neighborhood of the singularities, imposing a perturbation to the system 
up to the linear order, we obtain 
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3xy + 3—5- (5y + 3a;u<5i 

yd 



+ |'§_5T) + (33, 



35 / 2/ 

3x?/ + 3^- ] fa + 3yuSu 
x 6 I 



+ [ + ) + — zF- 61 ^ ( 34 ) 



t ( 3ua; + 3^-") fa + f-3tty - 3^ 



6^5h +6^« 2 , (35) 
a; ?/ 

4 \ / ;5 



3Z iy -3^ 5y + 6^ 



2 



y J j y 

3huSu, (36) 



W = (^ + 6 f)^+(-3y/ 2 -3ij^ 
3 ; 2 x + 3^fa-6^ 2 

x 6 J ar 

+ ihuSu, (37) 

where P is defined in (|26p. Before calculation of the eigenvalues of the linearized system around singularities, we 
must point out a key difference between the two singularities. The criterion for the stability of a small deviation must 
be with an argument which increases when the orbits go to the singularity. It is easy to check that the argument 
s = — ln(l + z) in (|2ip - (|2"5"l) is adapted to the attractor in late time universe. The dust attractor resides at the early 
universe, while the de Sitter attractor inhabits in late time universe. Thus, the linearized system about them should 
be described by arguments with opposite sign. If we use the same argument s, a positive definite eigenmatrix implies 
that the dust attractor is stable. 

For the dust attractor in the early universe, substitute the variables by the values given in (|2"5)) . (f2U)) . we obtain the 
eigenvalues of the linearized system 

3 3 

Ai =0, A 2 = 0, A 3 = 3, A 4 = -, A 5 = -. 

Hence, base on the above discussions, the dust attractor is quasi-stable. For the de Sitter attractor in late time 
universe, substitute the variables by the values given in (|30|) . (f3"Tjl . ([32]) . we obtain the eigenvalues of the linearized 
system 

3 

Ai = -6, A 2 = -6, A 3 =0, A 4 = 0, A 5 = --. 

which is also a quasi-stable attractor. 

Hence, eventually, the universe enters a de Sitter phase, which does not contain any singularities, which means the 
big rip disaster disappears naturally. The phantom behavior of the dark energy is only a transient phenomena. Also 
this result does not depend on parameter selection and initial condition because of attractor behavior. 
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The most significant parameter from the viewpoint of observations is the deceleration parameter q, which carries 
the total effects of cosmic fluids. Here q reads 



da 1 



l-8(^ + £) 



1 



Pi P2 Pi + P2 + Pn 



and density of Chaplygin gas u and density of dark matter v should satisfy 

x 2 (0)-y 2 (0) + u 2 (0) = 1. 



(38) 



(39) 



And then Friedmann equation ensures the spatial flatness in the whole history of the universe. Note that not all the 
examples in figure 1, 2 or 3 satisfies this constraint, since we do not introduce Friedmann equation there. Before 
analyzing the evolution of q with redshift, to obtain some asymptotic behaviors of the universe is useful. When 
z — > oo, q must go to 1/2 because both type I, II Chaplygin gas and matter behave as dust ; while when z — > —1, q 
is determined by 



lim q 



!im -[l-3(af 



(40) 



which agrees with the above analysis of the dynamical properties of this system. Here we carefully choose a new set 
of parameter which satisfies Friedmann constraint (|39p , and plot the deceleration parameter q in figure 2] We see for 
these three set of parameters the deceleration parameters are well consistent with observations. 



III. PERTURBATION GROWTH 



The Friedmann equation in the present model is given by f) 15[) . Now we consider the perturbation growth in this 
hybrid model. After the matter decoupling from radiation, for a region well inside a Hubble radius, the perturbation 
growth satisfies the following equation [80], 

6 + 2H5 - -^ Pm 5 = 0. (41) 



With the partition functions, 
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the perturbation equation (|4Tj) becomes, 



(In 5)" + (ln<5)' 2 + (^2 + (In 5)' = |f2 ro , 



where a prime denotes the derivative with respect to In a. f2 m and fide evolve as 

n m0 (i + z) 3 



a, 



y/A + B(l + zf -y/C + D{1 + zf 



(44) 

(45) 
(46) 



Ono(l + zf + ^A + B{l + zf - ^C + D{l + zf 
where denotes the present value of a quantity. The growth function defined in ([5]) is just (In a)'. Thus (|44|) generates 



f' + f 



1 3 

-(1 + fife) + -W de (Q m + - 1) 



/ — 2 ^ m ' 



(47) 



where we have used 



3 

— = - - [fi r , 



(1 + w)n de 



and 



Recalling ([5]), we derive the evolution equation for 7 from (|4T|) . 



3w de (l — fi m ) lnf2 



where we have used 

By using © and we obtain 

A 

Wde 



- W m - '^-"t + 3iu de f2 m (l/2 - 7) + 3 7 u> de - ~w de + 1/2 = 0, 



(lnO m )' = 3(1 - Cl m )wde. 



C 



_^A + B{l + zf y/C + D(l + z)« 

Expanding 7 about 1 — fl m , we get 

3 



^A + B{l + zf- ^C + D(l + zf 



- — + 2.4 x 1(T 2 (1 - fl m )(l - w de ){\ - Zw de /2)(1 - 6w de /5Y 

5 - Wde/(1 - Wde) 



(48) 
(49) 

(50) 
(51) 

(52) 
(53) 



Here we present some examples to illuminate the current numerical value of 7 in the hybrid Chaplygin gas model. We 
take the same examples as in figure 4. B = 0.3, D = 0.15, fi m = 0.28 with different A and C. For A — 5.92, C — 3.00, 
7 = 0.555; for A = 1.63, C = 0.300, 7 = 0.553; for A = 1.0JL C = 0.0300^7 = 0.550. The theoretical values of this 
model are well consistent with observations [HI, El, ES El El EE E3, El, 0, IzqI, [IE [zl III III • 



IV. CONCLUSION AND DISCUSSION 



Through careful analysis, we find a new kind of Chaplygin gas , called type II Chaplygin gas, whose energy density 
is negative. Then we present a hybrid Chaplygin gas model, in which type I Chaplygin gas and type II Chaplygin 
gas play the role of dark energy together. The EOS of dark energy crosses —1 naturally without introducing any 
gravity theories. 

In frame of standard general relativity and a spatially flat FRW universe, we study the dynamical properties of the 
present model. We find attractor solutions both in the early universe and the late time universe: The former greatly 
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mitigates the coincidence problem, while the latter overcomes big rip disaster. The stability about the singularities is 
also investigated. 

Some concrete examples of the deceleration parameter of this model are plotted. The result is well consistent with 
observations. Also because this model is definitely free of any extra difficulties from CMB decoupling, the structure 
formation can be studied in frame of it. 

We discuss the perturbation growth of this model. We find that the result is also consistent with observations. 
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